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Abstract 

For an electron gas with <5-function attraction we investigate the crossover from 
weak- to strong-coupling supercoductivity in two and three dimensions. We de- 
rive analytic expressions for the stiffness of phase fluctuations and set up effective 
XY-models which serve to determine nonperturbatively the temperature of phase 
decoherence where superconductivity breaks down. We find the transition tem- 
perature T c as a monotonous function of the coupling strength and carrier density 
both in two and three dimensions, and give analytic formulas for the merging of the 
temperature of phase decoherence with the temperature of pair formation in the 
weak-coupling limit. 
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1 Introduction 



The crossover from BCS superconductors to a Bose-Einstein condensate of tightly bound 
fermion pairs was first studied many years ago in Refs. in a model with 5-function 

attraction. This crossover has recently raised renewed interest [f|]- | 34fl , especially after 



the doping dependence of the critical temperature (T c ) and non-Fermi-liquid behavior 
above T c were combined in a phenomenological discussion of High-T c cuprates [pOj] , fl4} . 
In this paper we present a detailed study of the crossover based on a nonperturbative 
procedure in which derivative expansion is used to set up an effective XY-model whose 
well-known nonperturbative properties render information on the entire crossover regime 
of the above model. 

Physically, the most important distinctions between conventional (BCS) and strong- 
coupling (Bose-Einstein) regime lies in the fact that in the former only a small fraction of 
the conduction electrons is paired with the superfluid density involving all pairs, whereas 
in the latter practically all carriers are paired below a certain temperature T*, although not 
condensed. The temperature has to be lowered further below some critical temperature 
T c < T* to make these pairs condensed and establish phase coherence, which leads to 
superconductive behavior. 

In the model to be investigated in this paper, the crossover from BCS-type to Bose- 
type superconductivity will take place either by varying the coupling strength, or by 
decreasing the carrier density. 

Since in the BCS theory pair coupling is weak, it can well be described by mean-field 
theory for the pair fields. In the opposite limit of strong pair binding, on the other hand, 
superconductivity set in as a macroscopic occupationis of q = level and we are obliged to 
go beyond mean field level to describe such an effect. In three dimensions, crossover from 
BCS superconductivity to the Bose-Einstein condensation of tightly bound fermion pairs 
was first investigated in Ref. 0] by summing particle-particle ladder diagrams which corre- 
spond to Gaussian fluctuations around the mean field, in the functional integral formalism 
it was studied in Ref. |J. In both papers, fluctuation corrections were retained in the num- 
ber equation which was solved together with mean-field gap equation. In these papers, 
starting from fermionic system, gas of electron pairs was maped in the strong-coupling 
limit to the ideal Bose gas and critical temperature was shown to becomes asymptoticaly 
temperature of Bose-Einstein condensation of ideal Bose gas of particles with mass 2m 
and density n/2, where m and n are fermion mass and density. In this aproximation T c has 
an arificial maximum at intermediate coupling strength thus approaching limiting value 
in the strong-coupling limit from above. This artifact was removed in the generalized 
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self-consistent Green's function numerical approach in |J. 

In this paper we shall study the properties of collective modes with help of the lowest 
gradient terms governing the Gaussian fluctuations around the mean-field solution. These 
are most violent in the phase of the order parameter. Phase transitions in a system 
with these fluctuations are well understood in two and three dimensions from extensive 
studies of the XY-model. By setting up an equivalent XY-model we are therefore able 
to describe well the onset and disappearance of superconductivity in the entire crossover 
regime. In this way, we shall obtain simple formulas for the critical temperature T c as 
a monotonously increasing function of the coupling strength and carrier density in both 
two and three dimensions. In the weak-coupling limit, we give simple explicit formulas 
which show how the temperature of the XF-model transition converges to the transition 
temperature in the BCS theory. 

The weak- to strong-coupling superconductivity crossover in two dimensions was stud- 
ied via the Kosterlitz-Thouless theory in Refs. p7j, [p5fl . In Ref. p5] it was investigated 
within the same model as ours at a fixed carrier density, but only numerically^. The 
different properties of size and phase fluctuations was also exploited in Ref. |57| . 



2 Model 

The Hamiltonian of our model is the typical BCS Hamiltonian in D dimensions {% = 1) 
H = J2jd D xri(x)(-^-^ a (x)+9 /rf D o;V}(x)^(x)^(x)^(x), (1) 

where ^ CT (x) is the Fermi field operator, a =|, \ denotes the spin components, m is the 
effective fermionic mass, and g < the strength of an attractive potential g5(x — x'). 

The mean-field equations for the gap parameter A and the chemical potential \x are 
obtained in the standard way (see for example |j3~6|): 



- = — tanh^, (2) 



n = 77 = 77 V l-^Ltanhgp , (3) 




3 We shall see in Section 4 that these numerical results do not cover the entire crossover region, in 
particular, the merging of Tkt and T* in the weak-coupling region is missing - the effect that we show 
analyticaly in our paper. 
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where the sum runs over all wave vectors k, N is the total number of fermions, V the 
volume of the system, and 



k : 



2 



£ k = v /^ + A2 with ^=—- fi (4) 

are the energies of single-particle excitations. 

The 5-function potential produces an artificial divergence and requires regularization. 
A BCS superconductor possesses a natural cutoff supplied by the Debye frequency cup. 
For the crossover problem to be treated here this is no longer a useful quantity, since 
in the strong-coupling limit all fermions participate in the interaction, not only those in 
a thin shell of width up around the Fermi surface. To be applicable in this regime, we 
renormalize the gap equation in three dimensions with the help of the experimentally 
observable s-wave scattering length a s , for which the low-energy limit of the two-body 
scattering process gives an equally divergent expression |p|-|i~3|: 



m I 1 V— ft) 
4vra s g^V^V {) 



Eliminating g from (|5[) and (Q) we obtain a renormalized gap equation 



m 1 ^ 

4^~ 



— — tanh 
2E k 2T k 2 



(6) 



in which l/kpa s plays the role of a dimensionless coupling constant which monotonically 
increases from — oo to oo as the bare coupling constant \g\ runs from small (BCS limit) 
to large values (BE limit). This equation is to be solved simultaneously with @. These 
mean-field equations were first analyzed at a fixed carrier density in Refs. || and [[/J . Here 
we shall first reproduce some of the obtained estimates for T* and \i. 

In the BCS limit, the chemical potential fi does not differ much from the Fermi energy 
ep, whereas with increasing interaction strength, the distribution function broadens 
and fj, decreases. In the BE limit we have tightly bound pairs and nondegenerate fermions 
with a large negative chemical potential ^> T. In the strong coupling limit Eq. (|6]) 
provide us an estimation for T* - characteristic temperature of the thermal pair breaking 
0, whereas Eq. @ determines /i. From Eq . (|3|) we obtain that in the BE limit 
/i = —Eb/2, where E^ = l/ma 2 s is the binding energy of the bound pairs. In the BE 
limit, we can estimate that the pseudogap sets in at T* ~ Eb/2 log(£'fe/e j p) 3 / 2 . A simple 
"chemical" equilibrium estimate (/it = 2/i/) yields for the temperature of pair dissociation: 
Tkissoc — Eb/ \og(Eb/eF) 3 ^ 2 which shows at strong couplings T* is indeed related to pair 
formation ||, || (which in the strong-coupling regime lies above the temperature of the 
onset of phase coherence @-^T 
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The gap in the spectrum of single-particle excitations has a special feature [[| , |] , 0] 
when the chemical potential changes its sign. The sign change occurs at the minimum of 
the Bogoliubov quasiparticle energy where this energy defines the gap energy in the 
quasiparticle spectrum: 

£ gap = min(^ + A 2 ) 1/2 . (7) 

Thus, for positive chemical potential, the gap energy is given directly by the gap function 
A, whereas for negative chemical potential, it is larger than that: 



E 



A for // > 0, 

(/i 2 + A 2 ) 1 /2 f or /i<0 . 



In two dimensions, a nonzero bound state energy eo exists for any coupling strength. 
The cutoff can therefore be eliminated by subtracting from the two-dimensional zero- 
temperature gap equation |T0|-|T2| 



- = — V —, = = — / dz - 9 

9 2V k J + A 2 4vr J- xo ^TTz^ V ; 



where z = k 2 /2mA — xq, Xq = fi/A the bound-state equation 

1 1 1 m [°° d * (10) 

~~g~v\ k 2 /m + e ~ 2^ J-x Z 2z + e /A + 2x ' 

After performing the elementary integrals, one finds: 



^ = sfl +X 2 q -Xq. (11) 

In the next sectioins we will work at finite tempetrature, in doing so, we do not fix the car- 
rier density but assume the presence of a reservoir which provides us with a temperature- 
independent chemical potential // = ix{l/kpa s \T = 0). Such a fixed will be most 
convenient for deriving simple analytic results for the finite-temperature behavior of the 
system. In this fixed fi model carrier density becomes temperature dependent .0 For 



experimental measurements of \i see [38 



4 In Ref. [ p5[ , the temperature dependence of the chemical potential was calculated numerically for 
the entire crossover region within a "fixed carrier density model", where it turned out to be very small 
in comparison with the dependence on the coupling strength. 
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3 Phase fluctuations in Two Dimensions 
and Kosterlitz-Thouless Transition 



In this section we make use of derivative expansion to determine the relevant stiffness 
parameter for the study of phase fluctuations, which in two dimensions determines the 
temperature of the Kosterlitz-Thouless transition. In two-dimensional system, the phase 
fluctuations are most violent causing the strongest modifications of the mean-field prop- 



erties. The Coleman-Mermin-Wagner-Hohenberg theorem |39| forbids the existence of a 
strict long-range order, but there is quasi-long-range order manifesting itself in a power 
behavior of the correlation functions at all temperatures below Tkt- We shall neglect 
the coupling to the magnetic vector potential throughout the upcoming discussion, so 
that the phase coherence below T c can be of long or quasi-long range, unspoiled by the 
Meissner effect which would reduce the range to a finite penetration depth. 

The effective Hamiltonian from which we deduce the stiffness of the phase fluctuations 



was derived in [42], [25]. In this Section we summarize a few important aspects of it, with a 
reminder of its derivation given below. Writing the spacetime-dependent order parameter 
as A(x)e ld ( x \ where x denotes the four- vector x = (r, x) formed from imaginary time and 



position vector, the partition function may be written as a functional integral [[42], [40], [41 



Z(jjl,T) = J AVAVOexp T, A(x), d0(x))], 



(12) 



where 



T, A(x), d9(x)) = - [ dr f dxA 2 (x) - Tr logG -1 + Tr logG 1 (13) 
g Jo J 



is the one-loop effective action, containing the inverse Green function of the fermions in 
the collective pair field 



G~ 



-id T + r 3 +fA+ nA(r, x) 



T 3 



T 0(r,x) (V0(t,x)) : 



8m 



+ / 



2V 2 fl(r,x) zW(t,x)V 



4m 



+ 



2m 



(14) 



Here Ti, r 3 are the usual Pauli matrices, and Go = G\^a,9=o is added for regularization. 

Let us now assume that phase gradients are small. Then Q(fi,T, A(x),d8(x)) can be 
approximated as follows: 



ft(/z, A(x), d9{x)) ~ fi kin (/i, T, A, d9{x)) + fi pot (/i, T, A), 



(15) 
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with the "kinetic" term (see |42| 



OO 1 

n kin (/i,T,A,^(x)) =TTr^ -(</£)' 



n=l 



A=const 



(16) 



and the "potential" term 



fi pot (/i, T, A) = Q | d D o;A 2 - TTr log^T 1 + TTr logG ^ 



(17) 



A=const 



The latter coincides with the mean-field energy, determining the modulus of A(/j,,T) 
and thus the stiffness of phase fluctuations. The kinetic part i\m contains gradient 
terms whose size is determined by the modulus of A(/j,,T). Given the stiffness, one may 
immediately set up an equivalent XY-model. Both Q kin and Q pot are expressed in terms 
of the Green function of the fermions, which solves the equation 



and 



X(d6) = r 3 



id T e | (w) 2 



5772 



S(t,x) = 6(t)6(x) 
iV 2 9 «V0(r,x)V 



4m 



+ 



2m 



(19) 



The gradient expansion that we use to determine stiffness was first made in Ref. [42] at 
zero temperature. In Ref. ||25|| , the kinetic term £\m was calculated in two dimensions at 
finite temperature for arbitrary chemical potential retaining terms with n = 1,2 in the 
expansion fllpp. 
The result is 

T rP r 
^kin = — I dr I dPx 



2 Jo 



n(ji, T, A)id T 6 + J(/i, T, A(//, T))(W) 2 + T, A(/x, T))(«9 r #) 2 

(20) 

where J(/i, T, A) is the stiffness coefficient whose explicit form is 



1 



T r°° 



J (it, T, A) = - — n(a, T, A) , , 

4m 4tt 7-^/2T cosh 2 / x 2 + A 2 / 4T 2 



x + /i/2T 



(21) 



The other coefficents are: 

K(fi,T,A) 



m 
8^ 



1 



VJFTA 2 



tanh 



y//?+A^ 
2T 



(22) 



and n(/x, T, A) is the density of fermions @ which varies with temperature in our model. 
At the temperature T* where the modulus of A vanishes, also the stiffness disappears. 
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We are now ready to specify the effective XY-model governing the phase fluctuations. 
The model Hamiltonian corresponding to the gradient term in 9(x) is 



H= J -J dx[W(x)] 2 . (23) 

In contrast to the standard XY-model, the stiffness parameter is not a constant but 
depends on temperature via the solution of gap and number equations (0) and (|3]). Clearly, 
Kosterlitz-Thouless transition always take place below T*. For vortices of a high fugacity, 
the temperature of the phase transition is determined by the well-known formula 

Tkt = | J, (24) 

which follows from the divergence of the average square size of a vortex-antivortex pair. 
Since these attract each other by a Coulomb potential v(r) = 2it Jlog(r/r ), the average 
square distance is 

/•OO 1 

<r 2 >oc/ ^ r 2 e -^^)io g Wr„ )(X 1^ (25) 



'■() 



4-2ttJ/T' 



which diverges indeed at the temperature (|24|). In our case Tkt should be determined 
self-consistently: 

T KT = - J(n, T KT , A(/i, T KT )). (26) 

From pip, (0) and ( p4|) it is easily seen that Tkt indeed tends to zero when the pair 
attraction vanishes in which case A(T = 0) = 0. In general, the behavior of Tkt for 
strong and weak couplings is found by the following considerations. We observe that 
the particle number n does not vary appreciably in these limits with temperature in the 
range < T < T*, so that weak-coupling estimates for Tkt derived within the model 
with temperature- independent chemical potential (i.e., when the system is coupled to 
a large reservoir ) practically coincide with those derived from a fixed fermion density. 
Further it is immediately realized that in the weak-coupling limit A (Tkt, a0 /Tkt is a small 
parameter. At zero coupling, the stiffness J(/z, T K t, A(/i, T K t)) vanishes identically, such 
that an estimate of J at weak couplings requires calculating a lowest-order correction to 
the second term of eq.(pT|) proportional to A(T K t, aO/Tkt- Thus weak-coupling expression 
for stiffness reads: 

7C(3) A 2 (T) 
J(T) * W^-f^- (27) 
Equating this with the stiffness in (B3) we obtain the weak-coupling equation for T K t: 



e F A T 



■ KT 



1-^)- (28) 



4 V T* 
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where ep = {ji jm)n is the Fermi energy of free fermions. It is useful to introduce reduced 
dimensionless temperatures Tkt = Tkt/^f and T* = T*/eF which are small in the weak- 
coupling limit. Then we rewrite Eq. (^) as 



T, 



KT 



For small T* we may expand 



41 + 1/4T* 



T KT kT*- 4T*\ 



(29) 



(30) 



This equation shows explicitly how for decreasing coupling strength Tkt merges with T*. 
For weak coupling strengths, T KT behaves like 



T 



KT 



~ — A(0). 



7T 



(31) 



The merging of the two temperatures in the weak-coupling regime is displayed in Fig. [1]. 




10 20 30 40 50 

///A(O) 

Figure 1: Weak-coupling behavior of Tkt- The solid line is T*, the dashed line represent 
Tkt- 



Consider now the opposite limit of strong couplings. There we see from Eqs. (|26[) , 
(i) ©) an d (13) f° r ^kt, n(T,fi), and A(T,fi) that Tkt tends to a constant value. We 
can observe that in the strong-coupling limit A (Tkt) is always situated close to the zero- 
temperature value of A(T K t,/-0 ~ A(T = 0,/i). Taking this into the account we derive 
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an estimate for the second term in (0), thus obtaining the strong-coupling equation for 
Tkt' 



7T 



KT 



-n 



KT 



m 



71 



exp 



/i 2 + A 2 (T KT ,/i) 



KT 



(32) 



With the approximation A(Tkt,aO ~ A(T = 0,/i) we find that the first term in the 
exponent tends in the strong-coupling limit to a constant, A 2 (T K t, A i )/2/iT K T — ► —4, 
whereas the first term in the brackets tends to — oo, so that Eq. has the limiting form 



it n ( 1 
T K t - — 1 - - exp 
8ml 8 



(33) 



Thus for increasing coupling strength, the phase- decoherence temperature Tkt tends very 
quickly towards a constant: 



T, 



KT 



7T n 
8 m 



(34) 



In this limit we know from Eq. @ that the difference in the carrier density at zero 
temperature, n(T = 0), becomes equal to n(T = Tkt), so that our limiting result coincides 
with that obtained in the "fixed carrier density model": 



KT 



e F (n ) 



71 



5m 



(35) 



where we have inserted ep{n) = (7r/m)n for the Fermi energy of free fermions at the 
carrier density n Q = n(T = 0). 

From the above asymtotical formulas for weak- and strong-coupling limits we see 
that the temperature of the Kosterlitz-Thouless transition is a monotonous function of 
coupling strength and carrier density. The crossover takes place in a very narrow region 
/x/A(0) G (—1, 1). It is also observed in the behavior of three-dimensional condensation 
temperature T c of the gas of tightly bound, almost free composite bosons. In Refs. [0 , 
which include only quadratic fluctuations around the mean field (corresponding to ladder 
diagrams), T c was shown to tend to a constant free Bose gas value T c = [n/2C(3/2)] 2//3 7r/m, 
with no dependence on the internal structure of the boson. 

Here we find a similar result in two dimensions, where Tkt tends to a constant de- 
pending only on the mass 2m and the density n/2 of the pairs. No dependence on the 
coupling strength is left. The only difference with respect to the three-dimensional case 
is that here the transition temperature T c = Tkt is linear in the carrier density n, while 
growing like n 2//3 in three dimensions. Our limiting result ( pop agrees with Ref. |27[] and 
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25fl . There exists a corresponding equation for the temperature T* in the strong-coupling 
limit €o ^ ep: 



2 loge /e F ' 

4 Phase Fluctuations in Three Dimensions 



(36) 



In this section we discuss in a completely analogous way the fluctuations in three dimen- 
sions, where stiffness is, for small temperatures where A(T) is close to A(0): 

■fo(A*,T,A) = tU^T, A) - ^gi r^ —^+^ —, (37) 



4m ^' ' 7 16tt 2 T7- m cosh (VP + A 2 /2T) 
governing the phase fluctuations via an effective XY-model 

H = ^J d 3 x[W(x)] 2 . (38) 

The temperature of the phase transition in this model can reasonably be estimated using 
mean-field methods for the lattice 3D XY-Model ||: 

T& F ~ 3J 3D a, (39) 

a = 1/ n^ 3 is the lattice spacing of the theory || where is number of pairs. 

In the weak-coupling limit, the stiffness coefficent may be derived with the help of 
Gorkov's well-known method (setting T c « T*): 



7 ™ 3 A 2 

J 3D = T^C(3)^^, (40) 
48tt 4 ml* 1 



This is precisely the coefficient of the gradient term in the Ginzburg-Landau expansion, 
In the weak-coupling limit the two temperatures merge according to the formula: 

T c = T* — aT* 5/2 , (41) 

which contains a larger power of T* in the second term as well as a smaller prefactor 
a = (27r 2 ) 2 / 3 /2 w 3.65, as compared with the two-dimensional separation formula fl30|). 
The merging behavior is displayed in Fig. (|2]). 

In the strong-coupling limit of the theory where we have tightly bound composite 
bosons, the phase stiffness tends asymptotically to: 



= n_ _ 3V27rm 3/2 ex 
4m 167T 2 



V/^ 2 + A 2 



(42) 
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It obviously tends in this limit quickly to 

Tl 

Jbe = —- (43) 
Am 

An estimate for the critical temperature, obtained via the mean-field treatment of the 3D 
XY-model on the lattice reads in this limit: 
3 



T c = — 



2m 

This quickly tends from below to the value: 



^ 2/3 _J ^ T 3/2 3/2 f_VJFT^ S 

2 J n l/3 27/6^/2 J c m eX Pl 



T 

± c 



(44) 



^ DXV = l£ =« W H °- 2£ - <45) 
This result is very close to the temperature of the condensation of bosons of mass 2m and 
density n/2, which, as it was discussed in the introduction was obtained including the 
effect of Gaussian fluctuations into the mean-field equation for the particle number [§, ^ 
yielding^ 

T Bosons = [ n / 2 ((3/2)f 3 7T/m = 0.218e F . (46) 

5 In the study of the critical temperature via gaussian corrections to the number equation, mentioned 
in the introduction, the crossover behavior of the critical temperature shows an artificial maximum in 
the region of intermediate couplings J^, |5|, so that the limiting value is approached from above in the 
strong-coupling limit. This is in contrast to our three-dimensional Xy-model approach. 
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5 Conclusion 



We have studied the crossover from BCS-type to Bose-type superconductivity. For this 
purpose we have used the gradient expansion of the effective energy functional to set up 
an equivalent XY-model which allows us to investigate the onset of long-range order in the 
phase fluctuations. In two dimensions, we have given a simple analytic expression which 
shows how the resulting Kosterlitz-Thouless temperature Tkt at which quasi-long-range 
order sets in moves towards the pair-binding temperature T*, and merges with it in the 
weak-coupling limit. A similar expression was found in three dimensions. In the strong 
coupling limit we find that critical temperature tends in both two and three dimensions 
to a constant value as the chemical potential changes it sign. 

Let us finally remark that the separation of T* and T c has an analogy in the fer- 
roelectrics and magnets which also contain two separate characteristic temperatures, for 
example in the latter case — the Stoner- and the Curie-temperature. It also can be studied 
more precisely in a simple field theoretic model in 2+e dimensions with an O(n) symmetry 
for large n. In such a model, the existence of two small paramters e and 1/n has permitted 
us recently to prove the existence of two transitions, and to exhibit clearly their different 
physical origins [45|. 
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Note added in proof: 

Here and in our earlier preprint cond-mat/9804206 we criticize the discussion of the weak- 
coupling behavior of superconductive phase transition in two dimensions given in Ref. [25 ]. 
However, after our paper went to print in Phys. Rev. B 59, 12083 (1999)), the authors 
of [^| updated their discussion in their preprint |cond-mat / 9709034] version 2, which ap- 
peared in JETP 88, 685 (1999), and our criticism no longer applies. 
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